#Modelo Tabela 5 (coluna 1)
modelobeta_estado <- gamlss(PROP_EST ~ + NAC_PRESIDENTE * GENERO + NAC_SEC_GERAL * GENERO + NAC_TESOURARIA * GENERO + EST_PRESIDENCIA * GENERO + IDEOLOGIA + SITUACAO_REELEICAO + SITUACAO_REELEICAO * GENERO + EST_SEC_GERAL * GENERO + EST_TESOURARIA * GENERO + RACA, data = na.omit(BD_NOVO), family = BEZI())  
summary(modelobeta_estado)


#Modelo Tabela 5 (coluna 2)
modelobeta_estado_sem_out <- gamlss(PROP_EST ~  + NAC_PRESIDENTE * GENERO + NAC_SEC_GERAL * GENERO + NAC_TESOURARIA * GENERO + EST_PRESIDENCIA * GENERO + IDEOLOGIA + SITUACAO_REELEICAO + SITUACAO_REELEICAO * GENERO + EST_SEC_GERAL * GENERO + EST_TESOURARIA * GENERO + RACA, data = na.omit(BD_SEM_OUT), family = BEZI())  
summary(modelobeta_estado_sem_out)